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Prediction for transverse energy dependence of J/tp to 
Drell-Yan ratio in Au+Au collisions at RHIC energy was 
obtained in a model which assume 100% absorption of J/tp 
above a threshold density. The threshold density was ob- 
tained by fitting the NA50 data on J/^p suppression in Pb-fPb 
collisions at SPS energy. At RHIC energy, hard processes 
may be important. Prediction of J/^p suppression with and 
without hard processes were obtained. With hard processes 
included, J/ip's are strongly suppressed. 

PACS numbers: 25.75.-q, 25.75.Dw 



In relativistic heavy ion collisions J/^ suppression has 
been recognized as an important tool to identify the pos- 
sible phase transition to quark-gluon plasma. Because of 
the large mass of the charm quarks, cc pairs are produced 
on a short time scale. Their tight binding also make 
them immune to final state interactions. Their evolution 
probes the state of matter in the early stage of the col- 
lision. Matsui and Satz jl] predicted that in presence of 
quark-gluon plasma (QGP), binding of cc pairs into Jf-ip 
meson will be hindered, leading to the so called J/tp sup- 
pression in heavy ion collisions |^ . Over the years several 
groups have measured the J/4> yield in heavy ion colli- 
sions (for a review of the data and the interpretations see 
|,|). In brief, experimental data do show suppression. 
However this could be attributed to the conventional nu- 
clear absorption, also present in pA collisions. 

The latest data obtained by the NA50 collaboration ^] 
on J/-0 production in Pb-t-Pb collisions at 158 A GeV is 
the first indication of anomalous mechanism of charmo- 
nium suppression, which goes beyond the conventional 
suppression in nuclear environment. The ratio of 
yield to that of Drell-Yan pairs decreases faster with Et 
in the most central collisions than in the less central ones. 
It has been suggested that the resulting pattern can be 
understood in a deconfinement scenario in terms of suc- 
cessive melting of charmonium bound states [Q. 

In a recent paper Blaizot et al H] showed that the 
data can be understood as an effect of transverse energy 
fiuctuations in central heavy ion collisions. Introducing 
a factor e — ET/Erib) and assuming that the suppres- 
sion is 100% above a threshold density (a parameter in 
the model) and smearing the threshold density (at the 
expense of another parameter) best fit to the data was 
obtained. Capella et al |^ analysed the data in the co- 
mover approach. There also, the comover density has 



to be modified by the factor e. Introduction of this ad- 
hoc factor e can be justified in a model based on excited 
nucleons represented by strings Q . 

At a fixed impact parameter, the transverse energy as 
well as the number of NN collisions fiuctuate. The Fluc- 
tuations in the number of NN collisions were not taken 
into account in the calculations of Blaizot et al Q or in 
the calculations of Capella et al We have analysed 
[8| the NA50 data , extending the model of Blaizot ct al 
[5| to include the fluctuations in number of NN collisions 
at fixed impact parameter. The Et distribution was ob- 
tained in the geometric model, which includes these fiuc- 
tuations. It was shown that with a single parameter, 
the threshold density, above which the J/ip suppression 
is assumed to be 100%, good fit to the data can be ob- 
tained. In the present paper, we have applied the model 
to predict the Et dependence of J/tp to Drell-Yan ra- 
tio at RHIC energy. With RHIC being operational, it is 
hoped that the prediction will help to plan experiments 
to detect QGP. In the following, we will present in brief 
the model. The details can be found in 

In 1^ Et distribution for Pb-t-Pb collision was ob- 
tained in the Geometric model |^,|l^. In this model, 
Et distribution of AA collisions is written in terms of 
Et distribution in NN collisions. One also assume that 
the Gama distribution with parameters a and (3 describe 
the E'r-distribution in NN collisions. Et distribution for 
Pb+Pb could be fitted well with a = 3.46 ± 0.19 and 
f3 = 0.379 ±0.021. In fig. la, the experimental data along 
with the fit are shown. Parametric values of a and /3 
indicate that average Et produced in individual NN col- 
lisions is f3/a ~ .IGeV. This is to be contrasted with 
the average Et ~ IGeV produced in other AA collisions 

As mentioned in the beginning, we have assumed that 
above a threshold density ric, the charmonium suppres- 
sion is 100% effective Charmonium production cross- 
section at impact parameter b is written as, 

rf2^j/^/^2^^ ^j/V J d'^sTlff{s)T^"{s~h)S{h,s) (1) 
where T^'^^ is the effective nuclear thickness function. 



dzp{s,z)exp{-aabs / dz'p{s,z')) (2) 



with aabs as the cross-section for J/tp absorption by nu- 
cleons. The exponential factor is the nuclear absorption 
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survival probability, the probability for the cc pair to 
avoid nuclear absorption and form a J/V'- 5'(b,s) is the 
anomalous part of the suppression. Blaizot et al |^ as- 
sumed that J/ij] suppression is 100% effective above a 
threshold density (tIc), a parameter in the model. Ac- 
cordingly the anomalous suppression part was written 
as, 



S'(b, s) = e(nc - enp(b, s)) 



(3) 



where Up is the density of participant nucleons in the 
impact parameter space. 



np(b,s) =T^(s)[l 



-(TjvjvT'£!(b-s) 



[Ta 



(4) 



and £ — ET/Exib) — Ex/nfi/a is the modification factor 
which takes into account the transverse energy fluctua- 
tions at fixed impact parameter j^. This modification 
makes sense only when Up is assumed to be proportional 
to the energy density. Implicitly it was also assumed 
that the Et fluctuations are strongly correlated in differ- 
ent rapidity gaps. The assumption was essential as NA50 
collaboration measured Et in the 1.1-2.3 pseudorapidity 
window while the J/tJj^s were measured in the rapidity 
window 2.82 < y < 3.92 [Q. Strong correlation between 
Et fluctuations in different rapidity windows is explained 
in the Geometric model 

We calculate the J/V' production as a function of trans- 
verse energy, at an impact parameter b as, 



(5) 



in pA collisions point to a smaller value of dabs 4 to 
5 mb. The dashed line in fig. lb, corresponds to (Jabs—'^ 
mb and Uc = 3.42/m^^. However, it may be mentioned 
that (Jabs='^ uib does not allow a good fit to the pA and 
S-U data 
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FIG. 1. (a) Transverse energy distribution in Pb-fPb col- 
lisions, (b) J I'll) to Drell Yan ratio in Pb+Pb collisions as a 
function of transverse energy. 



where Pn{h^ Et) is the probability to obtain Et in n NN 
collisions, expression for which can be found in 

The Drell- Yan production was calculated similarly, re- 
placing charmonium cross section in eq. I by the Drell- 
Yan cross-section. 



dV^^'/d^fo ^ cT^^ I d^sTA{s)TB{s - b) 



(6) 



In fig. lb, we have compared the theoretical charmo- 
nium production cross-section with NA50 experimental 
data. The normalization factor a'^l'^ ja^^ was taken to 
be 53.5. The solid curve is obtained with Uabs=^A mb, 
and Uc = 3.8/m~^. Very good description of the data 
from 40 GeV onward is obtained. The model reproduced 
the 2nd drop around 100 GeV, (the knee of the Et dis- 
tribution). It may be noted that if the fluctuations in 
the NN collisions were neglected, equivalent description 
is obtained with threshold density ric — 3.75/to^, with 
smearing of the Q function at the expense of another pa- 
rameter. It is evident that in this model, the smearing 
is done by fluctuating NN collisions. Theoretical calcula- 
tions predict more suppressions below 40 GeV, a feature 
evident in other models also. It is possible to fit the 
entire Et range, reducing the J/i/'-nucleon absorption 
cross-scetion. Recent data ll^ on the J/^p cross section 



Present model can be used to predict Et dependence 
of J/ ip to Drell- Yan ratio at RHIC energy. Recent PHO- 
BOS experiment showed that for central collisions, 
total multiphcity is larger by 70% at RHIC than at SPS. 
We assume that Et is correspondingly increased [p^ . 
Accordingly, we rescale the Et distribution for Pb-|-Pb 
collisions and assume that it represent the experimen- 
tal Et distribution for Au-|-Au collisions at RHIC (small 
mass difference between Au and Pb is neglected). At 
RHIC energy the so-called hard component which is pro- 
portional to number of binary collisions appear. Model 
dependent calculations indicate that the hard compo- 
nent grows from 22% to 37% as the energy changes from 
yi=56 GeV to 130 GeV ||l||. However we choose to 
ignore the hard component in E'x-distribution. The mul- 
tiplicity distribution obtained by the PHOBOS collab- 
oration, in the rapidity range 3 <| rj |< 4.5 could be 
fitted well with or without this hard component. Indeed, 
it appears that the data are better fitted without the 
hard component ||l^. Global distribution e.g. multiplic- 
ity or transverse energy distributions are not sensitive 
to the hard component. In fig. 2a, filled circles repre- 
sent the "experimental" Et distribution for Au-f-Au col- 
lisions at RHIC, obtained by scaling the Et distribution 
in Pb-|-Pb collisions at SPS. The solid line is a fit to the 
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"experimental" £'T-distribution in the geometric model, 
obtained with a=3.09 and /3=0.495. Nucleon-nucleon in- 
elastic cross section (aNN) was assumed to be 41 mb at 
RHIC, instead of 32 mb at SPS Fitted values of a 
and P are interesting. With these values average Et pro- 
duced in individual NN collisions at RHIC is f3/ a— 0.16 
GeV. This is to be compared with the value 0.1 GeV 
for Pb-I-Pb coUisions at SPS. Average Et produced in 
individual NN collisions at RHIC is increased by 60%, 
compared to SPS energy. The apparent inconsistency is 
resolved if we remember that iTjvAr is increased by 30% 
from SPS to RHIC energy. 
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FIG. 2. (a)Expected Et distribution in Au-|-Au collisions 
at RHIC energy. The solid line is a fit in the geometric model 
(see text), (b) Expected Et dependence of J/ip to Drell-Yan 
ratio in Au-|-Au collisions at RHIC energy. The dotted (solid) 
lines are the predictions with (without) the hard scattering 
component, for two sets of parameters. The The dashed curve 
is the prediction of Blaizot et al 

While the Et distribution at RHIC may not be sensi- 
tive to the hard scattering component, the J/ip suppres- 
sion will be. In the model J/ijj suppression is 100% if the 
UpET / ET{b) exceeds the threshold density ric- With- 
out the hard component, transverse density in Au-|-Au 
collisions will be nearly same as that in Pb-|-Pb colli- 
sions at SPS. At RHIC, £'t(&) is increased, individual 
NN collisions produces more Et- Then anomalous sup- 
pression will set at larger Et (threshold density being 
same). However, with the hard component, the trans- 
verse density will be modified. For / fraction of hard 
scattering transverse density can be written as |Tq], 



.p(b, s) ^ (1 - /)np(b, s) + 2/n^'^'-'*(b, s) 



(7) 



with n^'^'^''(b,s) = aNNTA{s)TB{h-s). With hard 
component, transverse density is increased, as a result, 
anomalous suppression will set in at lower Et- 

In fig. 2b, we have presented the J/il) to Drell-Yan ra- 
tio obtained in the present model for Au-f Au collision 
at RHIC energy. We have presented the results for the 
two sets of parameters; (A) (Tafcs=6.4 mb, nc=3.8 fm~'^ 
and (B) <Tabs=4: mb, nc=3.44 fni"^- The solid lines are 
the prediction for J/iJj suppression neglecting the hard 
component in the transverse density, for the two sets of 
parameters (A) and (B) respectively. The dotted lines 
arc the prediction including the hard component. We 
find that without the hard component, J/V' suppression 
at RHIC is similar to that obtained at SPS energy. The 
2nd drop which occurred at 100 GeV in Pb-t-Pb colli- 
sions at SPS energy, now sets in around 180 GeV (knee of 
the Et distribution being around that energy). As men- 
tioned earlier, at RHIC, ET{b) is increased, and without 
the hard component the transverse density is nearly same 
as it was in Pb-|-PB collisions at SPS. Anomalous sup- 
pression then occurs at larger Et- As it was for Pb-|-Pb 
collisions at SPS, the two sets give nearly same suppres- 
sion for Et beyond 100 GeV. It is also an indication that 
anomalous suppression occur at larger Et- 

If the transverse density is modified to include the hard 
component, J/ip^s are strongly suppressed. In fig. 2b, the 
dotted lines presents the results obtained with 37% |l^ ] 
hard component. Nearly same suppression is obtained 
for set A and B. At knee, suppression is 6 times greater 
than corresponding suppression without the hard scat- 
tering component. Very strong suppression wash out the 
2nd drop, which was clearly visible at SPS, or in the pre- 
diction without the hard component. With hard compo- 
nent, as mentioned earlier, transverse density is increased 
and anomalous suppression sets in earlier. Blaizot et al 
jisf also predicted the Et dependence of the J/ip sup- 
pression. In fig. 2b, their result is shown (the dashed line). 
As mentioned earlier, fluctuations in number of NN colli- 
sions at fixed impact parameter was neglected. Also the 
model for Et production was not microscopic. However, 
at large Et, our prediction agrees closely with theirs. 

To summarize, prediction for J/ip suppression in 
Au-|-Au collisions at RHIC energy was obtained in a 
model which assume 100% absorption of J/tp above a 
threshold density. Transverse energy fluctuations as well 
as fluctuations in the number of NN collisions at flxed im- 
pact parameter were taken into account. The threshold 
density was obtained from the analysis of NA50 data on 
J/ip suppression in Pb4-Pb collisions at SPS energy. At 
RHIC energy hard processes are important. Prediction 
of -J/ip suppression, with and without the hard processes 
differ considerably. Without hard processes, predicted 
J/ip suppression at RHIC energy is similar to that ob- 
tained at SPS energy. The 2nd drop which occur at Et ^ 
100 GeV at SPS energy moves upward to 180 GeV. In- 
clusion of hard processes modifles the transverse density 
resulting in considerable larger suppression. Very large 
suppression washes out the 2nd drop, which was visible 
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at SPS energy or in the prediction without the hard pro- 
cesses. 
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